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A cornerstone of special relativity is Lorentz Invariance, the postulate that all observers measure 
exactly the same photon speeds independently on the photon energies. However, a hypothesized 
structure of spacetime may alter this conclusion at ultra-small length scales, a possibility allowed in 
many of the Quantum-Gravity (QG) formalisms currently investigated. A generalized uncertainty 
principle suggests that such effects might occur for photon energies approaching the Planck energy, 
Epianck = Mpi anc kC 2 ~ 1.22 x W 19 GeV. Even though all photons yet detected have energies 
E p h « E pianck, even a tiny variation in the speed of light, when accumulated over cosmological 
light-travel times, may be revealed by high temporal-resolution measurements of sharp features in 
Gamma-Ray Burst (GRB) lightcurves. Here we report the results of a study using the emission from 
GRB 090510 as detected by Fermi's LAT and GBM instruments, that set unprecedented limits on 
the dependence of the speed of photons on their energy. We find that the mass/energy scale for a 
linear in energy dispersion must be well above the Planck scale, something that renders any affected 
QG models highly implausible. 



I. INTRODUCTION 

Some Quantum-Gravity models (see P~H5] and ref- 
erences therein) postulate an inherent structure of 
spacetime (e.g. either a foamy, a discrete, or a lumpy 
spacetime) near the Planck scale Xpi = ^Gh/c 3 ~ 
10 _35 m. Since Lorentz symmetry is scale invariant 
(i.e. all scales are equivalent), the postulated existence 
of a special scale due to such QG effects can possibly 
lead to violations of Lorentz Invariance (LTV). One 
manifestation of such violations could be a dispersion 
in photon propagation, in which the speed of a photon 
in vacuo becomes dependent on its energy. 

A detection and measurement of LIV effects would 
be indisputably invaluable to our understanding of the 
nature of spacetime at extremely small scales. Fur- 
thermore, even setting an upper limit on the magni- 
tude of such effects can still prove very valuable, since 
measurements affected by physics at tiny Planck scales 
are rare and hard to perform, and can still be used 
to constrain and guide the relevant research. In this 
case, we have used high-quality measurements from 
GRB 090510 performed with both the GBM and LAT 
instruments on board the Fermi observatory to set 
the most stringent limits to date on the magnitude of 
such energy-dispersion effects. 

In the following, a brief introduction to the rele- 
vant formalism of LIV will be given, and some of the 
properties of GRB 090510 will be described. For more 
information on this study please refer to our Nature 
publication and its associated supplementary informa- 
tion [6]. 



Lorentz-Invariance Violation 

Consider two photons of energies Eh > Ei emitted 
simultaneously from a distant astrophysical source at 
redshift z. According to the postulated LIV effects 



these two photons will travel with different velocities 
and will arrive with a time delay At equal to: 
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Here Mqq is the "Quantum-Gravity mass", a pa- 
rameter that sets the energy scale at which such QG 
effects start to become important. Its value is as- 
sumed to be near the Planck Mass (Mpi = Uc/Xpi ~ 
10 19 GeV/c 2 ) and most likely smaller than it. The 
model-dependent parameter n sets the order of the 
LIV and is assumed to be one or two, corresponding 
to linear (At cx AE/M QG , with AE = E h -Ei~ E h ) 

and quadratic (At oc (Eh/MQo) 2 ) LIV respectively. 
The model-dependent parameter s n is equal to plus 
or minus one, and it sets the type of LIV: a posi- 
tive (negative) s n corresponds to a positive (negative) 
time delay or a speed retardation (acceleration) with 
an increasing photon energy. 

In this analysis we used the times and the ener- 
gies of the detected photons from GRB 090510 to set 
an upper limit on the strength of any LIV effects or 
equivalently a lower limit on Mqq. Since Mqq is ex- 
pected to be close and most likely smaller than the 
Planck Mass, a limit that is close to or, even better, 
over the Planck Mass is especially physically mean- 
ingful since it excludes most of the allowed parameter 
space, rendering any affected models highly implausi- 
ble. It should be noted that not all QG models that 
predict a spacetime structure actually require such an 
energy dispersion. Instead, many of the models are 
consistent with LIV but they do not explicitly require 
it. A model that actually requires LIV, and can be di- 
rectly constrained by such limits, is the stringy-foam 
model described in [7]. 
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GRB 090510 

Because of their short duration (typically with short 
substructure in the form of a series of narrow spikes) 
and cosmological distances, GRBs are well-suited for 
constraining LIV. In this study, we used measure- 
ments on the bright and short GRB 090510, which 
triggered both the LAT and GBM on May 10th 2009 
at 00:22:59.97 UT (hereafter all times are measured 
relative to this trigger time). Ground-based opti- 
cal follow-up spectroscopic data [8], taken 3.5 days 
later, exhibited prominent emission lines at a com- 
mon redshift of z = 0.903 ± 0.003. The GBM light 
curve (figures [Til] b,c; 8 keV - 40MeV) consisted of 
7 main pulses! After the first dim short spike near 
trigger-time (hereafter called the "precursor") , the flux 
went down to background levels. The main emission 
as detected by the GBM started at 0.53s and lasted 
for <0.5s. On the other hand , th e main emission as 
detected by the LAT (figures [ill] a, d-f; E>20MeV) 
started after the main GBM emission (0.65s vs 0.52s) 
and also extended to a significantly longer time scale 
(than the GBM emission) of about 200s. The emis- 
sion detected by the LAT extended to an energy of 
about 31GeV, which is the highest energy ever de- 
tected from a short GRB. The fact that this 31GeV 
photon was detected shortly after the beginning of 
the burst (~0.8s), and the fact that the LAT-detected 
emission exhibited a series of very narrow spikes (up to 
few tens of ms width) that extended to high (> tens 
of MeV) energies, allowed us to set stringent limits 
on the Quantum-Gravity mass. We used two inde- 
pendent methods to set these limits, described in the 
following section. 



II. LOWER LIMITS ON THE 
QUANTUM-GRAVITY MASS 

A. Using the arrival time of a single photon 

According to equation [l] if we knew the LIV- 
induced time delay (At) between two photons of ener- 
gies Eh > E[, then we should be able to make a mea- 
surement on Mqq (for some model dependent param- 
eters s n and n, and for a known measured z). How- 
ever, this time delay cannot be measured because we 
do not know the exact emission times and because 
we cannot safely assume that the two photons were 
emitted from the same location (hence assume that 
there were no extra propagation delays due to non 
co-location). However, what we can do is to first as- 
sume that the higher-energy (HE) photon was emit- 
ted some time during a lower-energy (LE) emission 
episode (which starts from time t s t ar t and extends up 
to at least the HE-photon detection time t HE ), then 
calculate a maximum time delay for the HE photon as 



At rnax = tnE — tstart, and finally calculate using equa- 
tion [l] a lower limit on the Quantum-Gravity mass. 

Since the end time of the lower-energy emission 
episode cannot be safely assumed, the method men- 
tioned above only works to constrain positive time de- 
lays (s n — +1). Also, this method is still valid even 
if the HE photon and the majority of the photons 
comprising the associated LE-emission episode are not 
emitted from the same exact location. Specifically, it 
only requires that there should be at least one LE pho- 
ton that was emitted from the same location as the 
HE photon and was detected during the lower-energy 
emission episode, something that is generally safe to 
assume. Lastly, this method is only very weakly sen- 
sitive to any possible spectral lags occurring intrinsi- 
cally at the GRB. Such lags have not been observed 
in short GRBs: they have been observed only in sub- 
McV energies (and not in the higher-energy range of 
interest here), and even when they are observed (i.e. 
in long GRBs), they are of the order of a typical spike 
width (few tens of ms), which is considerably smaller 
than the time delays typically used in this analysis. 
As it will be shown, spectral lags, if they are actually 
present in GRB 090510, would have a negligible effect 
on our main results and would certainly not affect the 
conclusion of this study. 

Setting stringent limits on the Quantum-Gravity 
mass requires an as-high-as-possible Eh [Eh usually 
^>E[ so Eh — Ei~Eh) and an-as small-as-possible At. 
For GRB 090510, we used the highest-energy photon 
detected (31GeV) for settings such limits. Even if 
another lower-energy photon corresponded to a more 
stringent limit, to be conservative, we only used the 
31GeV photon since its larger time delay makes it less 
sensitive to uncertainties in the choice of t star t and 
to any intrinsic spectral lags. Furthermore, its higher 
energy renders it less probable of being a background 
event than an alternative lower-energy candidate. In 
the following calculations, again to be conservative, we 
used values for the GRB's redshift and for the 31GeV 
photon's energy reduced by one standard deviation. 

The 31GeV photon had a reconstructed energy of 
30.5327;gyGeV (la confidence intervals), it was de- 
tected 0.829s post-trigger, and it coincided with one 
of the GBM pulses. Very thorough analyses confirmed 
this event to be a real photon (instead of a back- 
ground cosmic ray) associated with this GRB (e.g. 
from the diffuse galactic or extra-galactic emission). 
The properties of its associated event (absence of en- 
ergy deposition in the anti-coincidence detector, the 
signature of an electromagnetic shower in the instru- 
ment, etc.) and the results from our event classifica- 
tion algorithms strongly supported its gamma-ray na- 
ture. Furthermore, its strong temporal (detected dur- 
ing the prompt emission) and directional (less than 
one PSF from the reconstructed by Swift direction) 
coincidence with the GRB supported its association 
with this source. And independently of these consider- 



eConf C091122 



3 



2009 Fermi Symposium, Washington, D.C., Nov. 2-5 



3 



ations, simple analysis regarding the rate of expected 
background events (both cosmic- and gamma-rays) re- 
sulted in a probability of less than ~ 10~ 7 (> 5a) of 
the background producing such an event. 

The most conservative and of very high confidence 
assumption that can be made regarding the possi- 
ble emission time of the 31GeV photon was that it 
was not emitted before the beginning of the precur- 
sor (30ms before the trigger). For such an assump- 
tion, At = 0.829 + 0.03 = 0.859s and the associ- 
ated lower limit on the QG mass for linear energy 
dispersion (n — 1) is Mqg,x ^ 1.19xMpi an ck- An 
illustration that corresponds to this choice of t s t ar t is 
shown in sub-figure |HI^ i with the black solid (n=l) 
and dashed (n=2) lines. It should be noted that the 
validity of this assumption depends on the absence 
of a second undetected precursor before the detected 
one. Even though the Fermi LAT and GBM did not 
show evidence of any such second precursor, and even 
though no such two-precursor short GRBs have been 
detected, the existence of a second precursor start- 
ing at time t prej 2 before the start of the first detected 



precursor t 



pre,l 



-0.03s would decrease this most 



conservative limit by a factor of \ HE tprc - 2 . The pres- 
ence of any intrinsic spectral lags would correspond 
to a similar change in the limits. A negative spectral 
lag of ti a g ~ 0.01s (in which the 31 GeV photon was 
emitted ti ag time before the lower-energy associated 
events) would correspond to a decrease of our lower 
limit by a factor of At ^ t t ' ag ~ 1 since here ti ag « At. 

Our most conservative limit, described above, can 
be improved by considering that it is actually very 
likely that the 31GeV photon is associated with only 
the main emission starting at t star t = 0.53s instead of 
with the whole burst (i.e. assumming that the 31GeV 
photon is far more likely to have been emitted during 
the main emission episode instead someting before it) . 
This assumption is strongly supported by the fact that 
the 31GeV photon is expected to be emitted together 
with other lower-energy (sub MeV or MeV) photons 
that do not suffer from LIV delays and therefore mark 
its emission time. Using this assumption for the rele- 
vant emission time, the maximum time delay for the 
31GcV photon becomes considerably smaller and our 
limits significantly stronger MQG l i>3.42xMp/ onc fe. 
Similarly, we can make somewhat less conservative 
assumptions regarding the lower-energy emission in- 
terval and associate the 31GeV photon with the start 
of the >100MeV emission or the start of the >lGeV 
emission, yielding stronger limits, yet with less con- 
fidence. Such limits correspond to smaller time de- 
lays, therefore they are relatively more sensitive to 
uncertainties in the choice of t star t (because of now 
having to choose t s t a rt based on fewer photons) and 
to any possible intrinsic spectral lags. The results 
from all possible associations are shown in table I. As 
above, these limits are illustrated in sub-figure |III[i. 



It should be noted that all the results of this method 
yield strong limits on Mqq for linear LIV n = 1, while 
they give considerably less constraining limits for the 
quadratic case. Also, as mentioned before, the above 
results are for the case of speed retardation s n = +1 
since we cannot safely assume the latest emission time 
the 31 GeV photon was emitted. 

Finally, we note that the 31GeV photon arrives near 
the peak of a very bright and narrow spike in the soft 
gamma- ray lightcurve, which has a width of ~ 10-20 



ms (see figure III). Such an association is not secure 
since the 31GeV photon could have landed over the 
spike just by chance. However, just as an idea re- 
garding the magnitude of the limits that could be 
attainable provided such an association was secure, 
we calculate a significantly higher limit of Mqg,i ^ 
102x Mpi anc k for both speed retardation and accel- 
eration s„ = ±1 (since now we can constrain both 
the earliest and latest possible emission time of the 
31 GeV photon). Similarly, a weaker but indepen- 
dent and somewhat more robust limit on a possi- 
ble negative time delay (which constrains the super- 
luminal case, s n — —1) may be obtained from the 
~0.75GeV photon that is observed during the pre- 
cursor. This photon has a high probability of being 
from GRB 090510 (a chance probability of the back- 
ground producing such a photon of ^1.2xl0~ 6 cor- 
responding to ~ 4.6er), and the la confidence inter- 
val for its energy is 693. 6-854. 4MeV. For the case 
of speed acceleration (s„ = —1), we can calculate 
an maximum time delay of At<19ms and a limit of 
Mqg.i ^ 1-33 x Mpi anck (using again a reduced-by- 
one-standard-deviation photon energy in the calcula- 
tion). These two associations with an individual spike 
are shown in figure III i with the vertical dashed bands 
and in table I. 



B. Using the DisCan method 

An alternative method was also used for constrain- 
ing any linear in energy LIV effects. This method, 
called DisCan 9] (Dispersion Cancellation), extracts 
dispersion information from all the LAT-detected pho- 
tons (~30 MeV - ~30 GeV) and does not involve bin- 
ning in time or energy. It is based on the fact that 
any QG-induced time delays would be expected to 
smear the spiky structure of the lightcurve. The Dis- 
Can method applies different trial spectral lags (time 
delays that are inversely proportional to the energy) 
to the lightcurve until it finds the one that maximizes 
a measure of its "sharpness". The trial spectral lag 
that accomplishes this is equal and opposite in sign to 
the sum of any QG-induced and intrinsic-to-the-GRB 
spectral lags. 

The DisCan method was applied using the photons 
detected by the LAT during the interval 0.5-1. 5s post- 
trigger, the burst interval with the most intense emis- 
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Fermi/LAT: GRB 090510. T tat , = [8.30,9.50] s 
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FIG. 1: The value of the cost function (Shannon Informa- 
tion) versus the trial spectral lag value (8). The best value 
of 9, equal to zero, is annotated, and shown as a vertical 
solid line. The two dashed vertical lines left and right of 
the best value represent the 9 values which are 100 times 
less probable (-18 and +2 ms/GeV respectively) than the 
best 9 value, for the given data set. Thus the contained 
interval between the two dashed lines is an approximate 
error region, but does not reflect statistical uncertainties. 



sion. The results were not significantly sensitive to any 
variations on the stop times for the interval (±0.25s) 
or the energy upper limits (1, 3 & 100 GeV). The 
sharpness of the lightcurve was measured by using a 
cost function: the "Shannon Information" (eq. 11 of 
[9]). The value of this cost function is equal to the 
entropy (modulo a minus sign), and the difference of 
two values of the cost function are equal to the rel- 
ative probability that one value is more likely for a 
given data set. Figure [l] shows the value of this cost 
function versus the trial spectral lag value. The mini- 
mum value of the cost function (most probable value) 
is 0j^ 8 ms/GeV. The errors here correspond to the 
trial spectral lag values that are 100 times less prob- 
able than the best value of Oms/GeV, and are shown 
with the two vertical dashed lines in figure [l] 

However, this result does not comprise a definite 
measurement of the spectral lag value; instead it is 
just one of the results that are compatible with the in- 
herent uncertainties associated with our choice of time 
interval and energy range and with the limited statis- 
tics of the dataset. To estimate these uncertainties, a 
bootstrap analysis was performed, in which the Dis- 
Can method was applied to a randomized data set (a 
set produced by randomizing the association between 
the energies and times of the events). The reassign- 
ment destroyed any correlation with the energy, and 
therefore it removed any spectral lags. As a result, the 
application of this method to a randomized dataset 
should on average measure a spectral lag value that 
is equal to zero, while the width of the distribution of 
spectral lag values would be a measure of the statis- 
tical uncertainty of the measurement associated with 
our dataset. The results of the bootstrap analysis are 



FIG. 2: For the interval analyzed in Figure]!] to gauge un- 
certainty due to statistical variations we generated 100 re- 
alizations with the photon times randomized. This figure 
shows the distribution of the minimum trial spectral-lag 
value 9 m i n for these 100 realizations. From this distribu- 
tion, the 99%CL for 9 mm is 30ms/GeV. 

shown in figure [2] which shows the distribution of the 
measured spectral lags for each of the 100 random- 
ized datasets. The measured spectral lags lie within a 
value of < 30ms /GeV for 99% of the cases, and within 
a value of < 10ms /GeV for 90% of the cases. 

The combined result from minimizing the cost func- 
tion for the actual GRB 090510 dataset and from the 
bootstrap analysis on the randomized datasets is an 
upper limit on the energy dispersion equal to <30 
ms/GeV or Mq G > 1.22 x M Ptanck for linear energy 
dispersion of either sign s„ = ±1 at the 99% confi- 
dence level. This result is shown in table I, along with 
the results from the previous method. For reasons 
similar to those advanced in the previous subsection 
(improbability of intrinsic lags or fortuitous cancella- 
tion of quantum gravity and intrinsic dispersion) we 
take this result as an upper limit on LIV-induced dis- 
persion (i.e. we ignore any possible intrinsic spectral 
lags). 



III. CONCLUSION 

We have used high-quality measurements on GRB 
090510 performed by the GBM and LAT instruments 
on board the Fermi spacecraft to constrain tiny vari- 
ations on the speed of light in vacuo that are linear 
or quadratic to its energy. We used two independent 
methods to obtain conservative and unprecedented 
upper limits on the magnitude of such speed varia- 
tions. Our limits (Mqc,i >fewxMj> Imcfc ) strongly 
disfavour any models predicting such linear-in-energy 
variations in the speed of light. Limits of such strength 
(over the Planck mass) are especially physically mean- 
ingful, since they exclude almost all of the allowed pa- 
rameter space, rendering any affected models highly 
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t start 


Limit on \At\ 


Reasoning for t sta rt or method 


Ei 


Valid for s n 


Confidence 


Limit on Mqg,i 


Limit on Mqg,2 




(ms) 


(ms) 


used for setting the limits 


(MeV) 






{Mpianck) 


(10 10 GeU/c 2 ) 


fa) 


-30 


< 859 


start of any <MeV emission 


0.1 


+ 1 


very high 


> 1.19 


> 2.99 


(b) 


530 


< 299 


start of main<MeV emission 


0.1 


+ 1 


hie'h 


> 3.42 


>5.06 


(c) 


630 


< 199 


start of main>0.1 GeV emission 


100 


+ 1 


high 


> 5.12 


> 6.20 


(d) 


730 


< 99 


start of main>l GeV emission 


1000 


+ 1 


medium 


> 10.0 


> 8.79 


(e) 




< 10 


association with<l MeV spike 


0.1 


±1 


low 


> 102 


> 27.7 


(0 




< 19 


if 0.75GeV 7-ray from 1st spike 




-1 


low 


> 1.33 


> 0.54 


(g) 


\At/AE\ < 30ms/GeV 


Lag analysis of all LAT photons 




±1 


very high 


> 1.22 





TABLE I: Lower limits on the Quantum-Gravity mass scale associated with possible Lorentz Invariance Violation, that 
we can place from the lack of time delay in the arrival of high-energy photons relative to low energy photons, from our 
observations of GRB 090510. 



implausible. Any models predicting quadratic LIV are 
not significantly affected by our results. Our limits can 
be used to guide future research in Quantum Gravity 
and in general give insight on the nature of spacetime 
at minuscule Planck scales. 
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FIG. 3: Panel (a): energy vs. arrival time with respect to the GBM trigger time for the LAT photons that passed the 
transient event selection (red) and the photons that passed the onboard 7-ray filter (cyan). The solid and dashed curves 
are normalized to pass through the 31GeV photon and represent the relation between a photon's energy and arrival time 
for linear (n=l) and quadratic (n=2) LIV, respectively, assuming it is emitted at t s tart = —30ms (black; first small 
GBM pulse onset), 530ms (red; main < MeV emission onset), 648ms (green; > WOMeV emission onset), 730ms (blue; 
>GeV emission onset). Photons emitted at tstart would be located along such a line due to LIV time delays. Panels 
(b)-(f): GBM and LAT lightcurves, from lowest to highest energies. Panel (f) also overlays energy vs. arrival time for 
each photon, with the energy scale displayed on the right side. The dashed-dotted vertical lines show our 4 different 
possible choices for t start- The gray shaded regions indicate the arrival time of the 31GeV photon (on the right) and of a 
750MeV photon (during the first GBM pulse) (on the left), which can both constrain a negative time delay. Panels (b) 
and (c) show background subtracted lightcurves for the GBM detectors. Panels (d)-(f ) show, respectively, LAT events 
passing the onboard 7-ray filter, LAT transient class events with E > 100A/eV , and LAT transient class events with 
E > 1 GeV. In all lightcurves, the time-bin width is 10ms. 
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